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ABSTRACT 

We give a comprehensive statistical description of the Lya absorption from the 
intcrgalactic medium in a hydrodynamic simulation at redshifts 0.1 - 1.6, the range 
of redshifts covered by spectra of QSOs obtained with the Hubble Space Telescope. 
We use the ENZO code to make a simulation in a cube with side length 76 Mpc 
comoving usin g 75 kpc cells, fo r a Hubble constant of 71 km s _1 Mpc _1 . The best 
prior work, by Dave et all (|1999| ). used an SPH simulation in a 15.6 Mpc box with an 



effective resolution of 245 kpc and slightly different cosmological parameters. We use a 
popular cosmological model and astrophysical parameters that describe the ultraviolet 
background that phot oionizes and heats the gas. At redshifts z = 2 this simulation 
is different from data. iTvtler et alj ((20071 ) found that the simulated spectra at z = 2 
have too little power on large scales, Lya lines are too wide, there is a lack high column 
density lines, and there is a lack of pixels with low flux. Here we present statistics at 
z < 1.6, including the flux distribution, the mean flux, the effective opacity, and the 
power and correlation of the flux. We also give statistics of the Lya lines including 
the line width distribution, the column density distribution, the number of lines per 
unit equivalent width and redshift, and the correlation between the line width and 
column density. We find that the mean amount of absorption in the simulated spectra 
changes smoothly with redshift with DA(z) = 0.0102(1 + z) 2 - 2 52 . Both the trend an d 
absolute values are close to measurements of HST spectra by iKirkman et alj (|2007f ) . 
The column density distributio n and the line width dis tribution are also close to those 
measured from HST spectra bv lJanknecht et al. (2006;), except for the mode of the line 
width distribution which is smaller in the HST spectra, probably because of photon 
noise. Although some differences that we saw at z — 2 arc too subtle to be seen in 
existing HST spectra, overall, the simulation gives a good description of HST spectra 
at 0.1 < z < 1.6. 

Key words: quasars: absorption lines - cosmology: observations - intergalactic 
medium - numerical simulations. 



1 INTRODUCTION 

We are exploring the physical conditions in the IGM using 
hydrodynamic numerical simulations. Our goal is to make 
simulations that give simulated H I Lya absorption that is 
statistically identical to that measured in spectra of QSOs, 
where the best data give errors on the statistics of a few 
percent. We hope that such simulations will provide a com- 
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prehensive, accurate and reliable description of the physical 
conditions in the IGM. 

In previous papers we compared our opt ically thin 
ENZO simulations to data at redshifts z ~ 2. In I Jena et alj 
(|2005l ) (J05) we found that popular values of cosmological 
and astrophysical parameters gave simulated spectra that 
simultaneously match both the the mean amount of Lya ab- 
sorption and the distributi on of Lya line wid ths in ground 
based spectra of QSOs. In ITvtler et alj (|2007T ) however, we 
found that the same simulation does not simultaneously 
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match other statistics measured in data. The power spec- 
trum of the flux in the simulated spectra has significantly 
lower amplitude than we measure in data, by 50% on large 
scales and 20% on small scales. The lines in the simulated 
spectra are wider with larger 6-values, and the simulated 
spectra have different distributions for both the flux per 
pixel and the number of lines at various column densities 
than do the data. Here we use the same simulation to ex- 
plore redshifts 0.1 - 1.6 where the Lya line is visible in 
spe ctra obtained from space. 

iDave et al.l (|l999l ) used SPH simulations to explore the 
Lya forest at these low redshifts. One of their simulations 
used cosmological parameters close to those favoured today: 
Q, m = 0.4, fij, = 0.0473, fi A = 0.6, h = 0.65 and a$ = 
0.8. However, they used a box only 15.6 Mpc (we use a 
Hubble constant of 71 km s _1 Mpc _1 throughout the paper) 
on a side and an effective resol ution of 245 kpc, both far 
from ideal. It is now well known (|Penlll997l ; lBarkana fc Loebl 
12004 ISirkdl2007l ; iTvtler et alj|2007ft . that we need box sizes 
of the order of a hundred Mpc to contain the long mode 
density perturbations. In smaller boxes, we lack these long 
wavelength modes and we see significantly lower velocities 
and temperatur es at all densities in cluding those that make 
the Lya forest (|Tvtler et al.l [20071 ) . We also require small 
enough cells th at we can effectively reso lve the structures in 
the Lya forest. iMeiksin fc White] ((20041) f o und t hat the ideal 
cell size is < 40 kpc and in ITvtler et al.l (|2007l ) we confirm 
this. The 75 kpc cells that we use here are large enough 
that the line widths are significantly larger than they would 
be with much smaller cells. J05 showed that a cell size of 
150 kpc gives a poor representation of the IGM, sometimes 
showing trends with input parameters that are the opposite 
of what we see with small cells that resolve the Lya forest. 

In below, we briefly describe the simulation code and 
parameters we have adopted. In ^3]we describe the statistics 
of the flux in the Lya forest including the mean flux, flux dis- 
tribution and the flux power and autocorrelation function. 
In [J4] we describe the statistics of the lines in the Lya forest 
including the b- values, column densities and their correla- 
tions. In 3S]we give a brief comparison to data from HST. 
We discuss and conclude in [J6] 



2 ENZO IGM SIMULATIONS 

The numerical simulation that we describe in this 



paper (|Bodenheimer et al.l 120071 ) uses the Eulerian hy- 
drod yn amic cosmologic a l co d e ENZO (|Brvan et al.l 
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2007; 
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120071). The simulation contains both CDM 



and baryons in the form of gas, but no stars. The simulation 
was run with popular values for cosmological parameters: 
a flat geometry Qtotal = 1, comprising a vacuum energy 
density of Qa = 0.73, Q m — 0.27 (CDM plus baryons), a 
baryon density of fit, = 0.044, a Hubble constant of Ho = 
71 km s~ Mpc - and an initial power spectrum scalar 
slope of n 3 = 1.0 with a current amplitude of ag= 0.9. 

The ENZO code follows the evolution of the gas us- 
ing non-equil ibrium chemistry and cooling for hydro gen and 
helium ions (|Abel et al.lll997l ; lAnninos et ai1ll997l ). After 
reionization at 2 = 6, photoionization is provided using the 
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Figure 1. The photoionization rate from the UV emitt ed by 
QSOs and galaxies calculated by lHaardt fc Madaul feOOlh that 
we use in this simulation. The vertical axis is in units of 10~ 12 
photoionization per neutral H atom per second. 



lHaardt fc Madaul |2001) volume average UV background 
(UVB) from an evolving population of galaxies and QSOs. 
We show the rate of ionization of neutral hydrogen in Fig. [T] 

The simulations are optically thin so that all cells ex- 
perience the same UV intensity at a given time. We do not 
treat the transfer of radiation inside the volume, and we in- 
clude no feedback from individual stars, galaxies or QSOs 
except for that implied by the uniform UVB. 

As in J05, we use two astrophysical parameters to de- 
scribe the intensity of the UVB. The parameter 7912 is the 
rate of ionization per H I atom in units of the Haardt & 
Madau model discussed above, while X228 measures the heat 
input per He II ionization, again in units of the rate for the 
Haardt & Madau spectrum. In J05 we showed how we can 
use X228 t° adjust the widths of the Lya lines and the mean 
amount of absorption in simulated spectra, while 7912 will 
also adjust the mean amount of absorption. We can then 
choose values for X22S and 7912 that together match both the 
line widths and mean flux in data. The simulation that we 
discussed in this paper uses 7912 = 1.0 and X22 8 = 1-8 , and 
we called it simulation or box A in iJena et all (|2005D (J05) 
and ITvtler et al.l (|2007l ). 

W e initiate the simulations using an lEisenstein fc Hul 
\ 19991 ) power spectrum for the dark matter perturbations, 
that we insert at z = 60. The simulated volumes are all 
cubes with strictly periodic boundary conditions. Hence the 
power is input at a finite number of discrete wavenumbers. 

The amplitude of the power that we insert varies 
smoothly with wavenumber. We insert the amplitude ex- 
pected for the universe as a whole, with no random varia- 
tions associated with the finite box sizes. The power in the 
simulations is not adjusted to include the variations in mean 
density that we see in the universe on the scale of the boxes. 



Low z Lya forest simulations 3 



In this limited sense, the boxes contain information on scales 
much larger that their sizes. 

The simulation has a comoving box size of 76.8 Mpc 
and a comoving cell size of 53.25/i or 75 kpc comoving. 
The simulation has one CDM particle for each cell initially, 
and each dark matter particle has a mass of Mcdm — 9.5 x 
10 6 ft -1 M . 

Our simulation is barely large enough that the values 
we obtain at z — 0.1 will be reliable. The nonlinear scale can 
be quantified by calculating the wavenumber fcjvL such that 
(2^F ■fo NL P ZD {k)A-Kk 2 dk = 1 where P iD (k) is the three 
dimensional matter power spectrum. This indicates that the 
nonlinear scale is 14.11 Mpc (comoving) at z = 1.6 and 25.51 
Mpc (comoving) at z = 0.1. These correspond to 18.38% 
of the box size at z = 1.6 and 33.22% of the box size at 
2 = 0.1. In smaller boxes, the largest modes that fit inside 
the box have become non-linear before this low redshift and 
the results will be unreliable. 



2.1 Comparison of this simulation to data at z — 2 

In IJena et all (|2005T lf J05) we showed that at redshift 1.95 
various simulations using typical cosmological and astro- 
physical parameters gave a good match to both the mean 
flux transmitted in the Lya forest and the b-value distri- 
bution that we use to describes the line widths. We noted 
in Fig. 40 of J05 that the power spectrum of the flux for 
these simulations was similar to that of data, but we did 
not ex amine this in detail . 

In lTvtler et~ai] (|2007h we compared a wider variety of 
statistics at z = 2 to the simulations that we use in this 
paper, and to other simulations using the same input pa- 
rameters in smaller boxes. We found that if we matched the 
mean flux well, the line widths in the simulation were too 
big using as = 0.9 and X228= 1-8 times the standard UVB 
heating. If we instead use the smaller as values preferred by 
the WMAP 3 year data, the lines in the simulated spectra 
would be still larger. We now know that we could better 
match the line widths using X228 < 1-8, corresponding to 
a softer ionizing spectrum that gives less heating. However, 
from Fig. 39 of J05, we would simultaneously need to use 
7912 > 1 to preserve the match with the mean absorption, 
and hence it is unclear whether the combination of X228 
and 7912 parameters would be concordant with other data. 

We also found that the simulated spectra had too little 
power on both small and large scales. The small scales are 
an alternative way of noticing that the lines are too wide 
in the simulations. However, the larger scales are harder to 
adjust. We might use a larger as value or we might need to 
adjust the heating. 

We found that our simulations have too many lines with 
log Nhi< 14 cm -2 , a slight lack of lines with log Nhi= 14-15 
cm" 2 and a large lack with log Nhi> 17 cm -2 . We did not 
investigate whether these differences are themselves enough 
to account for the m issing power. We also found, confirming 
iBolton et alli|2007h . that our simulated spectra have too few 
pixels with a lot of absorption. 

We concluded that there were three classes of changes 
that might allow us to make simulations that exactly match 
data. 

First, there might be problems with the simulations. We 



use moderate resolution optically thin simulations that we 
hope will match the low density IGM. iTvtler et all (|2007f ) 
showed that using a much larger box does not help, while 
a factor of four smaller cell size made the difference larger. 
iKohler fc Gnedinl (|2007h showed that they obtain the cor- 
rect number of Lyman limit systems (LLS, with log Nhi> 
17.2 cm -2 ) at z = 4 using 2 kpc cells with radiative trans- 
fer. Perhaps a much larger version of their simulation might 
also match the entire column density distribution, the power, 
line widths and flux distribution. We should certainly add 
radiative transfer, especially that associated with the ion- 
ization of He II by discrete sources. This should give more 
realis tic temperatures a nd may improve the match to the 
data jBolton et all 2007 ). Early results such as lBolton et alj 
l|2004h and lPaschos et all (|2007l ) show that radiative transfer 
increases the temperature of the IGM at the cosmic mean 
density, the opposite of what we need. However radiative 
transfer might also give lower tem peratures at the h igher 
densities that make the Lya forest (|Bolton et al.|[2007r i. 

A second option is to change the astrophysics and espe- 
cially the heating of the IGM. A softer UVB will probably 
make spectra that are more like data. The third choice is to 
change the cosmological parameters, specifically we expect 
that a simulation with as> 0.9 would make spectra with 
more large scale power, more like the data. 

In summary, we know that our simulation is not ideal, 
but they are as good as any yet done, and a major improve- 
ment on those published at z < 2. 

When we compare simulations to data we would like 
to concentrate on the low density IGM that we can most 
readily simulate. We take this to mean log Nhi< 17 cm -2 
since the lack of lines becomes much larger at higher column 
densities and it is relatively easy to find Lyman continuum 
absorption in real spectra. At z = 2 the simulations contain 
some such lines, but factors of m any less than there are 
in real spectra (|Tvtler et all l2007t Fig. 15). Many of the 
statistics that we show below include lines of all column 
densities, including the few lines with log Nhi> 17 cm~ 2 in 
our simulated spectra. We will be analysing 2000 spectra 
at each redshift from 0.1 to 1.6. These spectra contain no 
lines with log Nhi> 17 cm' 2 at 2 = 0.1 and 26 such lines at 
2 = 1.6. 



3 STATISTICS OF THE FLUX IN THE Lya 
FOREST 

We make spectra of the flux transmitted through the IGM 
at Lya wavelengths using the code des cribed in J05 a nd fo l- 
lowing the improvements discussed in ITvtler et all (|2007h . 
We made 2000 spectra, each parallel to 2-axis of the box and 
beginning at a random location on the :ry-face of the box. 
We made 2000 spectra at each redshift from 0.1 - 1.6, in 
steps of 0.1. The random numbers that selects the x, y coor- 
dinates of each spectrum on the cube's face were different at 
each redshift and hence these spectra are approximately in- 
dependent, except that they all sample the same large scales 
structure in the box. 

We use a number of pixels per spectrum equal to the 
number of cells along a box side. We calculate the optical 
depth as a function of velocity in a spectrum, r(v), using 
equations in §4 of J05, and the flux from F = exp(-r), 
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Figure 3. The effective opacity and a power law fit. The fit is 
a function of the form tJ?* = A[l + z) a where z is the redshift, 
A = 0.0102 and a = 2.252. The curves above and below the solid 
line are best fits to the 2a upper and lower limits on the individual 
points. They do not show the range of possible fits. 

where F = 1.0 in the absence of absorption, and F — > at 
the centre of saturated absorption lines. The spectra from all 
simulations are smooth functions of v, with no discreteness 
from the cell size, because to determine r at each position 
we integrate over at least ±512 cells, to the left and to the 
right of the absorbing cell. Each spectrum is made at a fixed 
redshift. We convert from the Mpc per grid cell into velocity 
v in km s _1 of a spectrum using the H(z) for the chosen z, 
but we do not increment the z as we move down a sight 
line, and we do not change the H(z). Each spectrum is then 
frozen at the indicated redshift with no evolution. 

In Fig. [2] we show examples of the spectra. At lower 
redshifts we see systematically less absorption in total, fewer 
strong lines with large equivalents, and more regions with 
less than 5% absorption. 

3.1 Effective Opacity 

In Fig. [3] we show the mean effective opacity as a function 
of redshift. We calculate the effective opacity from the mean 
flux F using T e ff = —ln(F), as is normally done. The mean 
flux is the mean of all spectra that we made at each redshift. 
In Table Q] we give values for this and other statistics that 
have a single value at each redshift. 

We see that the effective opacity fall smoothly with de- 
creasing redshift, from 0.085 at z = 1.6 to 0.013 at z = 0.1. 
This is well fit by the power law that we show. 

3.2 Mean Amount of Absorption 

We define DA, the mean amount of absorption, as 
DA = 1 — F, where F is the mean flux in all the spectra that 
we made at that z. In Fig. |4]we show the DA as a function 



Table 1. Statistics of the Lya forest at redshifts 0.1 to 1.6. IT2 
is the ionization rate in units of 10 — 12 ionizations per H atom 
per second. r e j j is the mean effective optical depth, and DA the 
mean absorption. The values in the err column are 10 4 times the 
standard deviation of the DA values. Under a 2 ^ we give 10 3 times 
the variance in the mean flux in each spectrum. For <x 2 we list 
100 times the variance in the mean flux in each pixel in all the 
spectra F mo( ; e is the mode of the mean flux in each spectrum. 
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Figure 4. The mean absorption as a function of redshift. The 
solid line shows the fit through the points obtained from the op- 
tical depth fit Ffi t = exp(— -rf^) in Fig. [3] The upper and lower 
curves are fits to the 2a upper and lower limits on each point. 



of redshift. Since the effective optical depth is T e // << 1, 
we expect DA = 1 — F = 1 — exp(—T e ff) — T e ff as can be 
seen in the figure. The amount of absorption dec reases by a 
factor of 6.2 from 2 = 1.6 to 0.1. We know from I Jena et al.l 
l|2005h that had we used a simulation with a cell size of 
18 kpc instead of the 75 kpc used here, the mean flux would 
be reduced by about 0.02 at z = 2. 
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Figure 2. Spectra of the Lya absorption from the simulation. For each redshift we choose to show a spectrum with a mean flux that 
is approximately the same as the mean flux of all the spectra at that z. Although the velocity span of the box, and hence the spectra, 
decreases with the Hubble constant at low redshifts, we plot the same velocity range for all spectra. 



3.3 Distribution of Flux 

In Fig. [S] we show the distribution of flux per pixel for all 
spectra through the box, the flux pdf. The flux pdf is the 
number of pixels in each range of flux, flux bins, divided 
by the total number of pixels and then divided by the bin 
width, which is 0.05. We have 1024 pixels per spectrum and 
2000 spectra per redshift. Using this common definition, the 
values unfortunately depend on the bin width. Our bins are 
all the same width starting with the first bin at flux < F < 
0.05 with a centre at 0.025, and a last bin at 0.95 < F < 1.0. 
Our bin centres and the widths of two of ou r bins differ 
from lKim et al.1 (|2007r i and lBolton et all |2007l ) who choose 
a first bin at flux —0.025 < F < 0.025 with a centre at 
0, and a last bin at 0.975 < F < 1.025 with a centre at 
1.0. Since they divided all their bins by the same constant 
0.05 (Kim private communication), their flux pdf changes 
slope abruptly in the first and last bins, being a factor of 
2.0 below the value that we infer for flux zero and 1.0. We 
do not follow their definition because it is unexpected and 
hides trends. 

In Fig. [5] we see a smooth and continuous change in the 
shape of the flux PDf with redshift. As the redshift falls, the 
fraction of pixels with 0.95 < F < 1.0 increases while the 
fractions decline for all other flux value. At z = 1.6, 71% 
of pixels have F > 0.95 and only 1.2% have F < 0.05. At 
z = 0.1, 95.6% have F > 0.95 and 0.18% have F < 0.05. 
These values suggests that strong saturated lines that have 



F < 0.05 are 6.7 times less common at z = 0.1 compared to 
z = 1.6. This rate of change is similar to the factor of 6.5 
decrease in the effective optical depth and 6.2 in the mean 
amount of absorption. 

The fractional changes in the flux pdf values are similar 
for all redshifts and fluxes, except for F > 0.95 where the 
curves for the different redshifts cross. In detail, we see a 
faster relative change with redshift at the lower fluxes; the 
curves spread out more to the left of the plot. This sug- 
gests that strong deep lines that correspond to higher col- 
umn densities may evolve faster than do the weaker lines 
with smaller column densities, a trend that we will see more 
clearly in other statistics. 



3.4 Flux power 

We define the flux contrast as 

f(u) = (F/F) - 1, (1) 

where u is velocity and F is the mean flux from the all 
spectra that we made at a given z. We measure the one- 
dimensional flux power of each sight line in the z direction 
by computing the fourier transform of the flux contrast, and 
we present the average of the power from all sight lines at a 
given 2 in Fig. [6] The smallest k value that we show is that 
for one wave along an edge of the box. This is a constant in 
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Table 3. The flux power at different redshifts at selected values 
of the wavenumber k (s/km). The first column lists the redshifts 
while the next four columns list the value of the logarithm of the 
flux power at the wavenumber listed on top of each column. 
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Figure 5. The pdf of the flux per pixel. We use flux bins of width 
0.05 from 0.0 - 1.0, and we plot values at the bin centres, 0.025, 
0.075, 0.975. 



Table 2. The flux pdf at different redshifts. The first column is 
the redshift of the bin centre while the rest of the columns are the 
flux PDF values at different redshifts. All bins have a width of 
0.05, starting with the first bin from flux - 0.05. The pdf values 
in the bulk of the table are the fraction of pixels with flux value 
in the bin, divided by 0.05. The flux pdf values in each column 
sum to 1/0.05. 



flux redshift 





1.6 


1.3 


1.0 


0.7 


0.4 


0.1 


0.025 


0.236 


0.175 


0.110 


0.080 


0.057 


0.0354 


0.075 


0.057 


0.040 


0.024 


0.015 


0.011 


0.0069 


0.125 


0.054 


0.036 


0.023 


0.014 


0.009 


0.0062 


0.175 


0.054 


0.038 


0.023 


0.014 


0.010 


0.0063 


0.225 


0.057 


0.039 


0.025 


0.016 


0.010 


0.0070 


0.275 


0.059 


0.040 


0.026 


0.017 


0.011 


0.0066 


0.325 


0.066 


0.045 


0.030 


0.019 


0.012 


0.0072 


0.375 


0.073 


0.050 


0.034 


0.022 


0.014 


0.0086 


0.425 


0.084 


0.055 


0.040 


0.025 


0.015 


0.0091 


0.475 


0.095 


0.066 


0.047 


0.031 


0.019 


0.0100 


0.525 


0.110 


0.077 


0.056 


0.036 


0.023 


0.0128 


0.575 


0.131 


0.093 


0.066 


0.046 


0.030 


0.0164 


0.625 


0.162 


0.114 


0.081 


0.054 


0.037 


0.0194 


0.675 


0.207 


0.143 


0.105 


0.070 


0.047 


0.0261 


0.725 


0.269 


0.188 


0.140 


0.094 


0.063 


0.0348 


0.775 


0.376 


0.263 


0.193 


0.131 


0.091 


0.0523 


0.825 


0.563 


0.396 


0.291 


0.198 


0.140 


0.0849 


0.875 


0.972 


0.693 


0.511 


0.339 


0.250 


0.1543 


0.925 


2.174 


1.620 


1.197 


0.838 


0.598 


0.3848 


0.975 


14.201 


15.828 


16.970 


17.937 


18.551 


19.1107 



comoving Mpc, but in velocity units it falls with the decrease 
in the Hubble constant at low redshifts. 

We could compute the autocorrelation of the flux for 
a given velocity lag 8u directly from a flux spectrum as 
£f{8u) =< (F(u) - F)(F(u + 5u) - F) >, where F is the 
mean flux for all spectra at that z and the brackets refer 



to an average across the pixels of the spectrum. Instead we 
choose to obtain the autocorrelation from 

^f(Su) = ^ J p(ky kSu dk (2) 

where P(k) is the power of the flux contrast and we can use 
the line of sight averaged power directly because we have 
divided by the mean flux of the spectra at that z. In Fig. [7] 
we show the average autocorrelation profiles at each velocity 
lag from all lines of sight. 

In Figs. [6] and [7] we again see clear systematic trends 
with z. As the redshift decreases the power drops on large 
scales log k < —1.0 s/km, and increases on small scales. The 
power changes by less than a few percent at log k — —1.0 
s/km (a wavelength of 62 km s _1 ) for 0.1 < z < 1.6. The 
decrease in power on larger scales at lower redshifts comes 
from the decrease in the mean amount of absorption as the 
ionization of the IGM increases. The increase in the power 
on small scale (large k) at low z implies that the Lya lines 
are narrower at low z, which we will see is correct in H4.2I 
As we saw in J05, had we used a cell size of about 18 kpc, 
the flux power at z — 2 would fall by about 15% at k = 0.03 
and about 22% at k = 0.01 s/km. 

The autocorrelation falls smoothly with both redshift 
and distance. It does not show the crossing of curves that 
we saw for the power. 

3.5 Distribution of the mean flux per spectrum 

In Fig. [8] we show the distribution of mean flux per line 
of sight. We calculate the mean flux in each spectrum, and 
we plot the fraction of the 2000 spectra with a given mean 
flux. This statistic depends on the length of the spectra, 
which is 76 Mpc comoving. We would see less variation in 
the mean flux in longer spectra. As the redshift decreases the 
distributions shift to the right, to larger mean fluxes, and 
the distributions become narrower, with less variation from 
spectrum to spectrum. None of the spectra at z — 1.6 has a 
mean flux > 0.98, where as these are the most common mean 
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Figure 6. The mean ID flux power from all 2000 spectra at each 
redshift. We list the values in Table [3] 



Figure 8. The fraction of spectra, one per line of sight (LOS), 
with a given mean flux Fl. 
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Figure 7. The flux autocorrelation function against the velocity 
lag. We use arbitrary curves to connect the points at a given 
redshift. 



fluxes at z = 0.1. These distributions, especially the tail 
containing the occasional spectrum with a lot of absorption, 
reflect both the approximately exponential distribution of 
Lya line equivalent widths and the clustering of matter on 
large scales. 



4 LYMANa LINE STATISTICS 

W e obtain line s t atisti cs by fitt ing Voigt pro fi les as described 
in IZhang et all (|l997T ). As in iTvtler et all (|2004T l and J05 
(§5.1, 6.2) for many functions and statistics we consider only 
lines with 12.5 < log N H i< 14.5 cm -2 . 

We shall examine two sample of lines; a full sample of all 
lines with central optical depths r D > 1CP° and a sub-sample 
with t > 0.05. A sample that was largely complete to cen- 
tral optical depths r > 0.05 would require high enough 
resolution (FWHM < 10 km s _1 ) to completely resolve all 
the narrowest lines that would otherwise appear too shallow. 
The SNR per pixel must exceed 20 to give lcr in the centres 
of the limiting lines, barely enough to detect a line that is 
several pixels wide. To detect most lines with r = 0.05 we 
would prefer SNR > 60. The sample with r > 10~ 5 includes 
a large number of additional shallow very low column den- 
sity "lines" that are minor fluctuations in the flux. Existing 
HST spectra have far too low SNR to show such features. 



4.1 Line widths as a function of column density 

In Fig. [9] we show the b- values that quantify line width and 
H I column densities for Lya lines in spectra at redshifts 1.6 
and 0.1. We show only a fraction of all lines in our spectra 
at z = 1.6 to avoid saturating the plot. We show a curve 
on the plots that marks lines w ith line centre optical depth 
t = 0.05. From (|Spitzerlll97g| ), 



Nhi — \fiiTob/a \ , 



(3) 



where A D is the rest frame wavelength and a is the cross- 
section. The curve has equation 



b (km s" 1 ) = 15.0(0.05/r o )(A/ H //T0 12 cm- 2 ), 
and appears concave on the log-linear plots. 



(4) 
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Table 4. The distribution of H I column density values in Lya 
lines. We list values of log f(Nni) as a function of log Njji value 
in the first column and redshift given in the top row. We include 
all lines with line centre optical depth r a > 10 — 5 . We use column 
density bins of width 0.2 for log N m = 12.0-16.0 cm" 2 , and 0.5 
for log Njji= 16.0-18.0 cm -2 . We list N/A when there are no 
lines in the column density bin in any of our 2000 spectra at that 



loP" Nut 






redshift 






cm - ^ 


1.6 


1.3 


1.0 


0.7 


0.4 


0.1 


12.0 


-9.51 


-9.53 


-9.56 


-9.62 


-9.67 


-9.78 


12.2 


-9.74 


-9.79 


-9.85 


-9.92 


-9.99 


-10.15 


12.4 


-10.01 


-10.07 


-10.15 


-10.23 


-10.28 


-10.37 


12.6 


-10.29 


-10.38 


-10.47 


-10.54 


-10.61 


-10.71 


12.8 


-10.61 


-10.71 


-10.80 


-10.90 


-10.94 


-11.02 


13.0 


-10.92 


-11.04 


-11.11 


-11.21 


-11.26 


-11.36 


13.2 


-11.26 


-11.37 


-11.46 


-11.55 


-11.62 


-11.77 


13.4 


-11.61 


-11.71 


-11.81 


-11.88 


-11.95 


-12.16 


13.6 


-11.95 


-12.06 


-12.12 


-12.24 


-12.35 


-12.51 


13.8 


-12.32 


-12.43 


-12.52 


-12.65 


-12.80 


-12.92 


14.0 


-12.69 


-12.79 


-12.92 


-13.04 


-13.17 


-13.30 


14.2 


-13.08 


-13.17 


-13.36 


-13.51 


-13.63 


-13.58 


14.4 


-13.44 


-13.60 


-13.70 


-13.84 


-13.95 


-14.05 


14.6 


-13.81 


-13.97 


-14.14 


-14.14 


-14.33 


-14.33 


14.8 


-14.23 


-14.30 


-14.42 


-14.65 


-14.79 


-14.70 


15.0 


-14.61 


-14.66 


-14.83 


-14.97 


-14.88 


-15.02 


15.2 


-14.96 


-15.08 


-15.21 


-15.20 


-15.26 


-15.33 


15.4 


-15.39 


-15.35 


-15.56 


-15.50 


-15.51 


-15.60 


15.6 


-15.55 


-15.55 


-15.94 


-15.81 


-15.87 


-16.37 


15.8 


-16.10 


-15.93 


-16.26 


-16.23 


-16.16 


-16.20 


16.0 


-16.37 


-16.27 


-16.18 


-16.78 


-16.83 


-16.94 


16.5 


-17.19 


-17.27 


-17.26 


-17.40 


-17.38 


-17.69 


17.0 


-17.74 


-18.17 


-17.97 


-17.85 


-17.88 


N/A 


17.5 


-18.47 


-18.54 


-18.68 


-19.30 


-19.23 


N/A 


18.0 


-19.54 


-19.95 


-19.88 


-19.80 


N/A 


-19.66 



We also show on the plots a straight line that represents 
the increase in the minimum 6- value as a function of log Nhi- 
A visual examination of the plot shows no major changes in 
the distribution in the b-N plane with redshift, but we will 
see below that there are definite differences in the both the 
b and N distributions. 



4.2 Column Density Distribution 

In Table [4] we list values for the H I column density distri- 
bution, f(N), lines per unit z, per spectrum and p er unit 
linea r column density. With this standard definition (|Tvtlerl 
1 19871 ) the integral of f(N) over a range of linear (not log) 
column densities gives the number of lines in that column 
density range per unit redshift: 



dn(z) 
dz 



= / f(N)dN. 



(5) 



In Fig.QIJ]we show the f(N) distribution. A power law is 
a fair approximation to the distribution at all the N values, 
although the plot scale obscures significant deviations. We 
also show on the plot the sub-sample of the lines with r D > 
0.05. This sub-sample contains all lines at all except the 
lowest columns log Nhi< 12.5 cm -2 . On the plot of f(N) 
for z = 0.1 we show the fit to the f(N) at z = 1.6 in addition 
to a fit to the f(N) at z = 0.1. We see that there are fewer 
lines per unit z at lower z at all N values, and especially 



at high N values. The slope of the f(N) becomes slightly 
steeper at lower redshifts. As redshift drops, the number of 
lines with large column densities falls faster than the number 
with small columns. 



4.3 The doppler parameter distribution 

In Fig. 1111 we show distribution of lines wi dths for redshifts 
1.2 an d 1.6, together with the Hui-Rutlidge l|Hui fc Rutledgel 
Il999t ) function, 



dn b\ b% 
-=B-e X p{--) 



(6) 



where dn/db is the number of lines per km s _1 per unit 
2. The Hui-Rutledge function has only one parameter, b a , 
that describes the typical line widt h, where the peak o f the 
distribution is b peak = 0.94576 CT . In lTvtler et all (|2004 Fig. 
18) we showed that this function gave an excellent fit to the 
distribution of 6-values at z = 1.95. The function is also an 
excellent representation of the data at all low redshifts. In 
detail, the data have a wider distribution than the function, 
with more more lines at both small and large 6-values. 

We show in Fig. [T^] the b a parameter as a function of 
redshift. As previously, we show two different samples of 
lines; all lines with r > 10 -5 and the sub-sample with r > 
0.05. In both cases, the b a falls smoothly with redshift from 
~ 26.5 km s _1 at z = 1.6 to ~ 24.5 km s _1 at z = 0.1. 
The two samples give identical 6 CT values because there are 
very few lines with r < 0.05 and log Nhi> 12.5 cm" 2 , the 
mini mum that w e consi der when we evaluate /(6). We know 
from IJena et al.l (120051 ) that a simulation with a cell size 
of 18 km s _1 would give 6-values smaller by approximately 
1.0%. 

In Fig. [13] we show the /(&) for redshifts 0.1 and 1.6, 
each divided by the number of lines per unit redshift with 
column densities log Nhi= 12.5-14.5 cm -2 , as listed in listed 
in Table [5] The distributions have very similar shapes, al- 
though there are clearly more lines with smaller 6-value at 
low redshift. 

We list in Table [5] the b a and the median 6-value b me d 
(50% of values are smaller) as a function of redshift. We 
prefer the 6 CT with high SNR spectra because this statistic 
isolates the most common line width and it is insensitive 
to the precise shape of the 6-value distribution at low and 
high values. However in low SNR spectra, such as those from 
HST, the 6 met j may be advantageous because it is even less 
sensitive to the detailed shape of the distribution. In low 
SNR data, we expect a false excess of lines with low 6-values 
and a lack of lines with large 6-values. 



4.4 The equivalent width distribution 

Fig. [14] shows the distribution of the rest frame equivalent 
widths Wo for all lines at two redshifts. We list values in 
Table H 

Fig. H5lshows for two redshifts the number of lines with 
Wo > Wo, some minimum equivalent width. We find that a 
pair of straight lines can represent the points at each red- 
shift, hence the distributions are approximately broken ex- 
ponentials; ln(dn/dz) oc W\. 

In Fig. [16] we show how the number of lines with W > 
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Figure 9. Each point shows the 6- value and log Njji value (cm -2 ) of a Lya line in one of the spectra at z = 1.6 (left) and z = 0.1 (right). 
We use 2000 randomly sampled spectra at z=1.6 and 1749 spectra at z=0.1. The relative number of points between the two redshifts 
corresponds to the ratio between the logarithm of the number of lines at z=1.6 (85817) and z=0.1 (20637) in the column density range 
of the plot, log Nhi= 11.5 — 15.5 cm -2 . The dashed curve shows lines with central optical depth t = 0.05 from equation Q. Lines to 
the left of this curve have t < 0.05 and will be extremely hard to detect in HST spectra. The dashed straight lines is a linear fit to the 
minimum 6-valucs, b rn i n = (5.63 ±0.14)«o3(AT tf/ /10 13 ) + (16.99 ±1.78) km s" 1 This is a fit to the minimum in each of 10 bins of width 
0.2 from log Njji= 12.5 — 14.5 cm -2 at z=1.6. We show the same straight line at z=0.1 




Figure 10. The column density distribution f(N) at z — 1.6 and z = 0.1. The plus symbols show all lines with r > 10 5 , while 
the open squares are lines with r > 0.05. We show the same solid line on both plots. It is the power law, f(Nffj) oc N HI that 
best fits the number of lines with t > 10 — 5 and log Nhi= 12.5 — 14.5 cm" 2 in bins of 0.1, all at z = 1.6. The fit is log / (Nhi) = 
(-1.74 ± 0.04)Zog(iV//_r/10 12 ) - (9.28 ± 0.52). We extend the line to all log N H i values. The dashed line is the power law fit at z = 0.1, 
logf(N HI ) = (-1.84 ± 0.02)log(N HI /10 12 ) - (9.61 ± 0.34), for the same column density range. 
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Figure 11. The distribution of the b-values f(b) = d 2 n/dzdb of Lyo lines with logNjjj = 12.5 — 14.5 cm -2 . The histogram shows the 
values in 1 km s~ 1 wide bins while the curve shows the Hui-Rutlcdgc function. The area under the histograms and curves is the number 
lines per unit redshift with 12.5 < log Nhi < 14.5 cm -2 . 
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Figure 12. The redshift evolution of b a from Table \E\ from the 
simulated spectra. The solid line is a linear fit in the form = 
(1. 67 ± 0.12) *z + (23.09 ± 0.11) km s" 1 . The curves above and 
below the solid line are best fits to the 2a upper and lower limits 
on the individual points. They do not show the range of possible 
fits. 



Wo depends on redshift. We see that power laws are fair 
representations for all z and Wo values. 

In Fig. [17] we show the equivalent widths as a function 
of column density. For low column densities, lines are un- 
saturated and on the linear portion of the curve of growth, 



Figure 13. The doppler parameter distribution at two different 
redshifts. The distribution functions both have unit area. 

where W oc log Nhi, there is no dependence on the b- value. 
These are the most common lines with log Nhi ~ 14 cm~ 2 
and Wo ~ 0.3 A. However, for larger log Nhi values, the 
lines become saturated and the W depends on the fe-value. 
We see that when W > 0.5 A we have some lines with 
very large column densities. The simulation has factors of 
several too few lines with log Nhi> 17 cm~ 2 at z = 2, 
and hence it also lack lines with W > 0.5 A. For satu- 
rated lines with a fixed log Nhi value, we obtain a larger 
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Table 5. Line widths as a function of redshift in simulated spec- 
tra. We list the Hui-Rutledgc b a parameter and the median b- 
valuc b me£ j from lines with r > 0.05. The last column contains 
the log of the number of lines per unit z with log Nhi= 12.5- 
14.5 cm -2 . 

z b a b med \og(dn/dz) 



1.6 
1.5 
1.4 
1.3 
1.2 
1.1 
1.0 
0.9 
0.8 
0.7 
0.6 
0.5 
0.1 
0.3 
0.2 
0.1 



25.49 
25.64 
25.24 
25.32 
25.45 
25.13 
24.84 
24.62 
24.79 
24.12 
23.67 
23.86 
23.71 
23.35 
23.59 
23.38 



28.38 
28.25 
28.16 
28.13 
28.03 
27.97 
27.82 
27.36 
27.33 
27.08 
27.04 
26.90 
26.75 
26.72 
26.65 
26.62 



2.53 
2.50 
2.46 
2.43 
2.40 
2.37 
2.34 
2.31 
2.29 
2.25 
2.23 
2.21 
2.19 
2.15 
2.11 
2.08 



Figure 14. The equivalent width distribution, denned as 
f(W ) = d 2 n/dzdW , shown at z = 1.6 and z = 0.1. W is 
the rest frame equivalent width of an absorption line in units of 
A . The distribution is computed in bins of size 0.05 in equiva- 
lent width and is shown for W 0.025, the left edge of the first 
bin which is cantered at W = 0.05. The solid and dashed lines 
refer to z = 1.6 and z = 0.1 respectively. Compared to data at 
z = 2, our simulation lacks Lya lines with log Nhi > 17 cm -2 or 
Wo > 0.5 A. 
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Figure 15. The number of lines per unit z with equivalent width 
exceeding (Angstroms). We show the results for z = 1.6 
(above) and z = 0.1 (below). The straight lines with slopes s 
join at \Vq = 0.5. Compared to data at z = 2 our simulation 
lacks Lya lines with W > 0.5 A. 



Table 6. The equivalent width distribution, f(W a ), with W in 
A, at different redshifts. A value of N/A means there were no 
lines at that equivalent width for that redshift. 



Wo 



redshift 





1.6 


1.3 


1.0 


0.7 


0,1 


0.1 


0.05 


3.490 


3.380 


3.299 


3.199 


3.160 


3.064 


0.10 


2.923 


2.826 


2.738 


2.655 


2.583 


2.403 


0.15 


2.613 


2.513 


2.428 


2.330 


2.208 


2.087 


0.20 


2.419 


2.289 


2.225 


2.112 


1.976 


1.801 


0.25 


2.220 


2.137 


1.997 


1.882 


1.787 


1.602 


0.30 


2.094 


2.003 


1.926 


1.727 


1.531 


1.585 


0.35 


1.997 


1.873 


1.662 


1.625 


1.501 


1.492 


0.40 


1.858 


1.804 


1.643 


1.523 


1.328 


1.301 


0.45 


1.807 


1.696 


1.582 


1.464 


1.372 


1.301 


0.50 


1.719 


1.591 


1.464 


1.376 


1.420 


1.253 


0.55 


1.616 


1.542 


1.302 


1.256 


1.201 


1.000 


0.60 


1.435 


1.325 


1.258 


0.981 


0.957 


0.801 


0.65 


1.263 


1.154 


1.041 


0.926 


0.860 


0.567 


0.70 


1.004 


0.957 


0.852 


0.584 


0.355 


0.199 


0.75 


0.889 


0.717 


0.740 


0.583 


0.355 


0.500 


0.80 


0.733 


0.341 


0.464 


0.487 


0.258 


-0.278 


0.85 


0.517 


0.216 


N/A 


0.062 


-0.042 


N/A 


0.90 


0.371 


0.392 


-0.189 


-0.115 


-0.343 


-0.278 


0.95 


0.216 


0.040 


-0.013 


-0.416 


-0.343 


0.023 


1.00 


-0.151 


0.040 


N/A 


N/A 


N/A 


-0.278 



Wo value when the b-value is larger which allows more of 
the atoms to absorb photons because they are more widely 
distributed in velocity. Hence the lines with W ~ 1 A and 
log Nhi~ 14.5 cm" 2 have large 6-values while those with 
Wo ~ 0.5 A and log Nhi> 16 c m" 2 have small b-values. We 
know from lMisawa et all (|2004 Figs. 3 & 5) that at z = 2.5 
some such lines have 6- values smaller than those of any lines 
with log Nhi~ 14 cm" 2 . Note that the total W from all lines 
in the simulated spectra with log Nhi> 17 cm" 2 is negligi- 
ble and hence these lines have a negligible effect on most 
statistics. We do not know how many lines we might see in 
a simulation with much smaller cell size and self-shielding. 
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slope = 1 .40±0.20 



z = 



0.2 0.3 
log(1 +z) 



Figure 16. The redshift evolution of the number of lines ex- 
ceeding equivalent width W o =0.1, 0.3 or 0.7 Angstroms. We also 
show power laws dn/dz oc (1 + z) 7 with slopes 7. 



4.5 Number of Lines per unit redshift 

We show the number of lines per unit redshift in Fig. 1 181 for 
both the large sample of lines with r > 10~ 5 and the sub- 
sample with t > 0.05. Note the dramatic difference! The 
large sample shows a factor of a few more lines at lower red- 
shifts, while the sub-sample, that is a better representation 
of data, shows the opposite. The number of li nes is clearly 
very sensitive to the definition of the sample (|Zhang et al.l 
1 19971 ). yet in both cases the change with redshift is very well 
described by a single power law. 

In Fig.[l9]we examine the redshift evolution of the num- 
ber of lines with very high column densities, log Nhi> 16 
and log Nhi> 17 cm" 2 . The lines with log Nhi> 17 cm" 2 
evolve more quickly. We should remind the reader that our 
simulation has factors of many too few lines with these col- 
umn densities because it is optically thin and the cell size is 
factors of many too large. 

We show in Fig. [201 the number of Lya lines above a 
cutoff column density at redshifts of 1.6 and 0.1. We see 
an excess of lines above the power law fit at 2 = 1.6 for 
log Nhi> 16 cm -2 . This excess is less at lower redshifts and 
it almost disappears by z — 0.1. 



5 COMPARISON TO DATA 

Our simulations have known limitations. They lack high col- 
umn density lines compared to simulations that have much 
smaller cell size and allow self-shielding. We also know that 
at z — 2 there are differences between simulated and rea l 
spectra that show up in most statistics l|Tvtler et al.ll2007t ). 

When we compare to data we would ideally match the 
properties of the observed spectra in detail, including the z 
range, signal-to-noise ratio, spectral resolution and contin- 
uum errors. In this section, we compare statistics of the Lya 
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Figure 17. The equivalent widths (Angstroms) and H I column 
densities (cm -2 ) of Lye? lines at 2 = 1.6 and z = 0.1. We show 
all lines in the plot range from all 2000 simulated spectra at each 
redshift, except for four with log Njji> 18.2 cm -2 at z = 1.6. 



forest in real and simulated spectra. We will ignore spectral 
resolution and the signal-to-noise ratio of the real spectra, 
although we will attempt to match the redshift and column 
den sity ranges. 

IKirkman et~all (|2007h measure the mean amount of ab- 
sorption in the Lya forest at low redshifts using HST spec- 
tra of 74 QSOs. In Fig. [21] we show t hat our DA values are 
very close to the observed ones from IKirkman et al.l (|2007l . 
Table 5), where we do not consider the dip in DA at z = 1.4 
to b e significant. 

Ijanknecht et alj (|2006l ) fit the column densities and b- 
values to lines in HST spectra of 9 QSOs at z em < 1.94. 
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Figure 18. The number of lines with central optical depth 
t > 10~ 5 (plus symbols) or r > 0.05 (open squares). We also 
show power laws dn/dz = A(l + z)t , with A = 3.439 ± 0.004 and 
7 = -0.820 ± 0.014 for t d > 10" 5 (above) and A = 2.24 ± 0.01 
and 7 = 1.06 ± 0.03 for r D > 0.05 (below). 



Figure 20. The number of lines per unit redshift as a function 
of the minimum column density at z = 1.6 and z = 0.1. We 
show both the large sample of lines (plus symbols) and the sub- 
sample with t > 0.05 (open squares). The solid straight line is 
a power law fit at z = 1.6 (plus symbols) to log Njji* = 12.5 — 
14.5 cm -2 (vertical dotted lines). The dashed straight line is the 
corresponding fit to the z = 0.1 points. 
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Figure 19. The redshift evolution of number of lines with 
log Njji> 16 cm -2 (open circles) and log Njji> 17 cm~ 2 (filled 
circles). The straight lines are power laws dn/dz cx (1 + z)" 1 . 
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They fit the column density distribution with a power law 
/(iV) = ^log Nhi 13 and find no significant change of the 
power with redshift, just as we do. They find a power law of 
f3 ~ 1.57 ± 0.05 consistent with our value of ~ 1.62. 

The median b- values from the HST spectra, 28 km s~ , 
and from our simulation, 28.5 km s , are also very similar, 



Figure 21. Comparison of the mean am ount of absorption, DA 
values, in HST spectra (error bars) from iKirkman et alj i2007l) , 
and in the simulated spectra (crosses). We have not adjusted the 
the simulation output or spectra in any way. 
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Table 7. Comparis on of Lya lines in HST spectra 
( Janknccht et al.l l2006h and in our simulated spectra. The 
HST spectra have have 270 lines at z = 0.5 — 1.0 and 595 lines 
at z = 1.0 — 1.5. /3 is the index on the power law equation of the 
distribution of column densities, f(N). The median line width is 
t>med (km s^ 1 ) and b a is the Hui-Rutledge parameter. 
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P 




bmed 




bo 






obs 


sim 


obs 


sim 


obs 


sim 


0.5 - 1.0 


1.57 ±0.05 


1.62 


29 


28.3 


18.9 


25 


1.0 - 1.5 


1.58 ±0.04 


1.63 


28 


28.5 


23.5 


25.5 



see Table [7] However, our b a values are larger than those in 
the HST spectra. We have two comments on the b a values. 
It is well known that samples of lines from spectra with sig- 
nificant noise, including all HST spectra, include an excess 
of lines with low 6-values. Common lines with the lowest 
column densities enter the sample when photon noise gives 
them extra depth in on e or two pixels: seelKirkman fc Tvtlerl 
(| 19971 . Fig. 10 & 11) or lMisawa et"al1 (|2004 Fig. 3). We saw 
in J05 that the simulation also makes lines that are too wide 
at z = 2 and that smaller cells would reduce the line widths 
but not enough t o resolve the differe nce. We discussed this 
issue at length in lTvtler et all (|2007l) . 



6 CONCLUSIONS 

Some statistical descriptions of the Lya forest vary greatly 
over the redshift range z = 0.1 — 1.6 while others change 
little. Statistics that change greatly derive from the mean 
amount of absorption: the flux pdf, the power and auto cor- 
relation of the flux, and the amplitude of f(N) and the 
number of lines per unit z. The statistics that change lit- 
tle relate to the line widths: the shape (not the amplitude) 
of the b- value distribution, /(&), the b a value and the b-N 
distribution. 

The changes in the Lya forest with redshift are all ex- 
ceedingly smooth, with a constant rate of change across the 
redshift range 0.1 - 1.6. There are no special redshifts and 
no need for "broken" power laws or changes in slope. We can 
readily predict the values of parameters at larger or smaller 
redshifts. Specifically, we do not see the jump in the evolu- 
tion of the number of lines numbe r of lines per unit z, as 
reported by I Wevmann et al.l i|l998h . Rather, our results are 
very close to the constant rate of evolutio n found in recent 
analys is of larger samples o f HST spectra bv lJanknecht et all 
(|2006T ) and lKirkman et al.l (|2007T ). 

We have found good agreement between the main statis- 
tics of the Lya forest lines in HST spectra and in our simu- 
lated spectra. These statistics include the mean amount of 
absorption, the change in the absorption with redshift and 
the median fe-value. The b^ is smaller in HST spectra than 
in our simulations. This difference might be real, but given 
the agreement in the median b-value, it is more likely a re- 
sult of the low SNR of the HST spectra giving a false excess 
of narrow lines. 

The agreement between the Lya forest in our simula- 
tions and in HST spectra is better than we had anticipated. 
We already know that there are differences between the sim- 
ulation and data at z = 2 but these are mostly too subtle to 
be detected in the existing HST spectra. The HST spectra 



with high enough resolution to resolve Lya lines are few and 
they have very low SNR by ground based standards. While 
we had adjusted the astrophysical parameters of the simu- 
lation to approximately match data at z = 2 we made no 
such adjustments for any lower redshifts. The evolution of 
the Lya forest in the simulation is controlled by the change 
in the gas density, the growth of structure and especially the 
change in the ionisation and heating by the UVB. We know 
that the UVB intensity was adjusted by Haardt and Madau 
to match a variety of data at various redshifts, but they did 
not explicitly match the low redshift IGM using a simula- 
tion of the Lya forest. They clearly did an excellent job of 
calculating the UVB at low redshifts. In this limited sense, 
our simulation is predicting the Lya forest at low redshifts, 
and it is doing much better than we had expected. This 
implies that the simulation comes close to representing the 
main physical process in the IGM, those responsible for the 
appearance and evolution of the Lya forest. 
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